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Recent cosmological data analyses hint at the presence of an extra relativistic energy component 
in the early universe. This component is often parametrized as an excess of the effective neutrino 
number A'efi over the standard value of 3.046. The excess relativistic energy could be an indication 
for an extra (sterile) neutrino, but early dark energy and barotropic dark energy also contribute 
to the relativistic degrees of freedom. We examine the capabilities of current and future data to 
constrain and discriminate between these explanations, and to detect the early dark energy density 
associated with them. We find that while early dark energy does not alter the current constraints 
on A^eff, a dark radiation component, such as that provided by barotropic dark energy models, 
can substantially change current constraints on A^eff, bringing its value back to agreement with the 
theoretical prediction. Both dark energy models also have implications for the primordial mass 
fraction of Helium Yp and the scalar perturbation index n^. The ongoing Planck satellite mission 
will be able to further discriminate between sterile neutrinos and early dark energy. 



I. INTRODUCTION 

The precision of theoretical modelling and ob- 
servational measurements of the cosmic microwave 
background (CMB) temperature and polarization 
anisotropies from satellites and ground based exper- 
iments [1|-1J] has opened the exciting possibility of 
addressing key questions about the nature of dark 
energy, dark matter and primordial inflation. Interest- 
ingly, a first hint of some new physics may be showing 
up in the high redshift universe from CMB and Big 
Bang nucleosynthesis (BBN) data. Recent analyses (for 
example [l|, 043) may be hinting at the need for an 
extra, dark, relativistic energy component. 

If further data confirms this, it could suggest new dark 
matter such as a sterile neutrino [5[ or a decaying par- 
ticle [sl-fToj. or nonstandard thermal history Or it 
could indicate that dark energy does not fade away to the 
^ 10~^ fraction of the energy density at CMB recombi- 
nation that is predicted by the cosmological constant. 
Indeed, some proposed particle physics explanations for 
dark energy involve scaling fields [l^ [l3| with an early 
energy density which is a constant fraction of the energy 
density of the dominant component. Another possibil- 
ity is that the evidence for the extra relativistic compo- 
nent may be signaling the presence of a "dark radiation" 
component in the early universe, as predicted by certain 
higher dimension braneworld scenarios 



ponent. One defines the effective number of neutrinos 
iVeff through : 
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Any of these would be exciting extensions to the stan- 
dard, concordance model. Uncovering new degrees of 
freedom would be of great importance, and distinguishing 
between the possible origins could give valuable insight 
into physics and cosmology. The extra, dark contribution 
to the total relativistic energy density in the early Uni- 
verse is generally phrased in terms of the energy density 
of neutrinos - the known dark, (early) relativistic com- 
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where pi, is the neutrino energy density and is the 
CMB photon energy density, with value today p^^o 
4.8 X 10"'^'* gcm~'^. In the standard model, with three 
masslcss neutrinos with zero chemical potential and in 
the limit of instantaneous decoupling, A^cff = 3. The 
inclusion of entropy transfer between neutrinos and the 
thermal bath modifies this number to about A^off = 3.046 
at the CMB epoch (see e.g. ji]). 

The recent analysis of [J] that combined CMB data 
with measurements of baryon acoustic oscillations (BAO) 
and the Hubble constant reported an excess of the rel- 
ativistic neutrino number, A"eff = 4.6 ± 0.8 at 68% c.l., 
disfavoring the standard value at about two standard de- 
viations. This is compatible with previous analyses [HIB" 

a. 

Another possible high redshift discrepancy, also sensi- 
tive to relativistic degrees of freedom, involves primor- 
dial '*IIc measurements compared with the predictions of 
Big Bang nucleosynthesis (see, e.g., While stan- 

dard BBN, assuming a value of the baryon-photon ra- 
tio of 7? = (6.19 ± 0.15) X 10^1° as measured by CMB 
data [l|, predicts a primordial Helium mass fraction 
Yp = 0.2487 ± 0.0002, current observational measure- 
ments prefer a larger value of 3^, = 0.2561 ± 0.0108 [T6j 
andFp = 0.2565±0.0010(stat.)±0.0050(syst.) [H. Since 
the primordial "^He mass fraction is largely determined by 
the neutron to proton ratio at the start of BBN, Yp is sen- 
sitive to the value of the expansion rate and so through 
the Fricdmann equation to the overall energy density at 
temperature ~ 1 MeV , e.g. in relativistic particles. From 
one has approximately Al^ ~ 0.013(A"off — 3) for 
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l^cfF — 3| < 1. Thus, an increase in the effective number 
of relativistic degrees of freedom could remove the ten- 
sion between BBN and the measured "^He abundances. 

These new relativistic degrees of freedom (RDOF) 
could be due, for example, to a fourth (or fifth), ster- 
ile neutrino (a postulated neutrino species that does not 
participate in weak interactions). Such a hypothesis is 
worthwhile testing and may also be cornpatible with re- 
cent neutrino oscillation results (e.g. jl9l - l2l| ). This origin 
would have no direct relation to the question of cosmic 
acceleration and dark energy. 

However, the signal could also arise from the class of 
early dark energy models. Scalar fields from dilatons in 
field theory and moduli in string theory are generally pre- 
dicted to possess scaling properties, so that they would 
evolve as radiation in the radiation dominated era and 
contribute a constant fraction of energy density. Extend- 
ing probes of dark energy to high redshift is an impor- 
tant frontier and detection of its effects would provide an 
invaluable guide to the physics behind cosmic accelera- 
tion. Dark energy that is significantly present not only 
in the late universe but also at early times is called early 
dark energy (EDE; see, e.g., [H,!!^). Furthermore, some 
theories, typically involving higher dimensions, predict a 
"dark radiation" component in addition to a cosmologi- 
cal constant. 

All three origins increase RDOF but have different 
evolutions in the energy density as radiation domination 
wanes, and hence will have the expansion history and N^s 
differing as a function of time. In this paper we consider 
the effects of contributions to A^eff from both neutrino 
and dark energy components, individually and together, 
and analyse the constraints imposed by current and fu- 
ture cosmological data. The main motivations are 1) to 
investigate how the current indication for RDOF would 
translate into a signal for dark energy at high redshift, 
and 2) to examine how adopting a dark energy compo- 
nent that is not negligible at high redshift would impact 
the stability of future possible conclusions from Planck 
about RDOF. 

As a product of this analysis, we update the current 
constraints on EDE with recent data. Previous analy- 
sis have placed constraints on EDE using the available 
cosmological datasets and forecasting the discriminatory 
power of future CMB probes like Planck (see e.g. |24l - 
|28|). Here we revise these constraints by using updated 
datasets, by enlarging the parameter space through in- 
cluding shear viscosity in EDE perturbations, and by 
considering the possible degeneracies with sterile neutri- 
nos. Section [IT] describes in more detail the models used 
to give extra A'^^ff . In Sec.|lll]we present the results of our 
analysis for the different cases, and Sec. IIVI discusses the 
conclusions about the ability to constrain and distinguish 
the various physical origins. 



II. NEUTRINOS, EARLY DARK ENERGY, 
AND DARK RADIATION 

We consider three models that contribute to RDOF: 
sterile neutrinos, early dark energy, and a barotropic dark 
energy model that produces a dark radiation component 
in the early universe. Accounting for each possible con- 
tribution, the RDOF translated into an effective number 
of neutrinos is : 

= A^eff + AiV,!^^ + AA^fff , (2) 

where N^g is the number of neutrino species (including 
extra sterile neutrinos) , AN^^^ is the contribution com- 
ing from early dark energy, and AA^^^ is the contribution 
from barotropic dark energy. When the components do 
not behave completely relativistically, the effective num- 
bers may be functions of redshift; for EDE the contribu- 
tion is constant only well before matter-radiation equal- 
ity, while barotropic dark energy behaves as a relativistic 
component at all times during and before recombination. 
In the following sections we describe our modelling for 
these three components and then their physical signa- 
tures in CMB power spectra, before proceeding to the 
ability of cosmological data to constrain and discriminate 
among them. 



A. Relativistic Neutrinos 

For the purposes of exploring a deviation from the 
standard value of Aoff = 3.046, we first assume ther- 
mal (so the factors in Eq. ([T]) hold), massless sterile 
neutrinos which give a time-independent contribution to 
Acff according to Eq. ([T]). We indicate the total neu- 
trino contribution to Acff by N^^, which is not necessar- 
ily equal to A'cff any more. Current cosmological data 
bound the mass of extra (thermal) sterile neutrinos to be 
m^,s < 0.5 eV for Acff > 4 at 95% c.l. (0, 0). Such 
massive sterile neutrinos may also be compatible with 
recent neutrino oscillation results (e.g. |19l421| ). however 
considering them as massive has negligible impact on the 
constraints on TVeff (again, see [l^) and we treat these 
neutrinos as massless in what follows. Neutrinos with 
~keV masses, sometimes considered for sterile neutrinos, 
are non-relativistic at recombination and contribute lit- 
tle to A'cff at that time. Note that, while we refer N^^g to 
sterile neutrinos in what follows, other relativistic back- 
grounds (for example, gravitational waves) produce iden- 
tical effects on cosmology. See HO] for a decaying particle 
scenario. Any such model that lacks significant contri- 
bution to late time energy density is for our purposes 
equivalent to sterile neutrinos. 

The usual case, e.g. in is to analyze the constraints 
on N^g in the absence of an early dark energy density, in 
which case A"cff = -'XTs- consider in this article Nj^g 
in the presence of an early dark energy density. 
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B. Early Dark Energy 

Early dark energy is the name given to a dark compo- 
nent that in the recent universe acts to accelerate expan- 
sion, but which retains a non-negligible energy density 
at early times (e.g. around recombination, or earlier). To 
keep an appreciable energy density at early times, the 
equation of state would not be negative, but at or near 
that of the background equation of state. We adopt the 
commonly used form [22| 



= ^"^0 :L ^ +»e(l 



w{a) = 



(ilnfidc(a) 
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where fide (a) is the fractional energy density and w{a) 
the equation of state of EDE. The factor acq/ (a -I- acq) in 
Eq. ^ comes from ilr{a)/[l — 57dc(a)] where ilr is the 
fractional radiation energy density (specifically excluding 
any EDE) and Ooq is the scale factor at matter-radiation 
equality. Here 51^^, and are the current dark energy 
and matter density, respectively, and a flat Universe is 
assumed so fi^j + O"^ = 1. The present equation of 
state w{a = 1) = wq- The energy density r2de(a) goes 
to a finite constant fig in the past, in both the matter 
dominated and radiation dominated eras, indicating a 
scaling solution. 

The dark energy equation of state w{a) follows three 
distinct behaviours: w ~ 1/3 during radiation domina- 
tion, w sa during matter domination, and w ~ wq in 
recent epochs. An accurate fitting formula for the time 
variation of the EDE equation of state during the recent 
universe is Wa = 5^,^ [2j1, where w{a) = wq + Wail — a) 
fits observational quantities at z < 3 to 0.1% accuracy. 
We extend the model by modeling the EDE clustering 
properties through the effective sound speed = 6p/6p 
and a viscosity parameter c^,^ that describes the possi- 
ble presence of anisotropic stresses (see e.g. HI] and ref- 
erences therein). In what follows we assume these clus- 
tering parameters as constant and consider two cases: 
= (^vis = corresponding to a relativistic ori- 



gin, and cj — 1, 



as expected in the case of 



a quintessence scalar field. For simplicity we consider 
Wq = —1 since the low rcdshift data are consistent with 
a cosmological constant and viable EDE models have lit- 
tle time variation there. 

Regarding the contribution to A'off, the EDE scaling 
behavior indicates the energy density will behave as a rel- 
ativistic component until the epoch of matter-radiation 
equality and so AN^^^ will be constant by Eq. ([1]). 
However, since w then deviates from 1/3 toward as 
the EDE later behaves more non-relativistically, AN^^^ 
will grow. Translating the EDE density into an effective 
"neutrino" number AN^^^ through Eq. ((!]) yields : 
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This is clearly a redshift-dcpcndent quantity since as 
the EDE equation of state begins to evolve differently 
from radiation the density ratio will vary with time. In 
the limit a <^ a^q, 



(a < aeq) = 7.44 



(6) 



since Pde/p-, = (Pde/Prad) (Prad//07) and /9rad/P7 = 1.69 

for three neutrino species. 

In Figurefflwe plot AN^ff^{a), for fie = 0.05. As we 
can see, at early times EDE behaves like a RDOF com- 
ponent with a constant value of AN^^^ « 0.39. How- 
ever this value increases at later times, when EDE starts 
to mutate into a matter-scaling component, reaching 
AiV^^^ « 1.6 at recombination. This time dependence 
will be a crucial element in discriminating between EDE 
and a sterile neutrino contribution A^^ff . Having AA^cft in 
EDE models smaller at BBN than at recombination helps 
ease the discrepancy between the lower value expected 
(3.046) and that derived from CMB data. Furthermore, 
the larger value of AA'eff at recombination means that the 
constraints on EDE from CMB anisotropies will translate 
to tighter bounds on AN^ff^ at BBN, and hence on Yp, 
than those in the neutrino RDOF case. 




FIG 1" ^ atEde 

a, for — 0.05 (the results scale nearly linearly for smaller 

values). Note the strong time dependence near recombina- 
tion. 



C. Barotropic Dark Energy 

As a model with some characteristics of each of the 
previous two, we consider a case containing dark radia- 
tion, that is a component whose energy density evolves as 
p (X but does not interact electromagnetically. Such 
terms arise in higher dimensional theories with multi- 
ple branes that induce a Weyl tensor contribution to the 
energy- momentum tensor |14l |. Here we choose a more 
conventional model with interesting properties, arising 
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from the barotropic class [3l| where the pressure is an ex- 
phcit function of the energy density. Barotropic models 
were shown to be highly predictive, reducing the cosmo- 
logical constant fine tuning problem by rapidly evolving 
through an attractor mechanism toward w = —1 js^ . 
Their equation of state, and hence energy density evolu- 
tion, is wholly determined by their sound speed, through 

w' = dw/d\na = -3{l + w){cl-w). (7) 

This gives for any constant Cg a time dependent equation 
of state 

w = [c2Sa-3(i+-=) - l]/[Ba-3(i+-?) + 1] , (8) 

where B = (1 + wo)/{cl — wq). 

Being interested in relativistic degrees of freedom, we 
choose cl = 1/3 (and indeed > 1/3 would violate 
early radiation domination). This leads to a surprisingly 
simple solution: 

Pbaro(a) = Poo + Cprfia~'^ , (9) 

where = (3iJ^/87rG)(l - 17™ - OT^,o) and C = 
ilf I (1 — 17^). This acts like radiation at early times, with 
a constant fractional energy density fi^ during the radi- 
ation dominated era. At late times it has a constant ab- 
solute energy density poo- It basically looks like the sum 
of a cosmological constant and dark radiation, despite 
having no explicit cosmological constant. As expected, 
at early times = 1/3 and at late times w rapidly ap- 
proaches — 1. We take wq = —0.99 (since w — —1 is only 
reached asymptotically), and c^^^ ~ 1/3 to match the 
other cases. 

To clearly state the main practical differences of our 
three models: extra neutrino species give a constant 
contribution to A'eff and negligible contribution to late 
time energy density as well as no acceleration; standard 
early dark energy gives a time varying contribution to 
TVeff as well as late time energy density and acceleration; 
barotropic dark energy has the third interesting combina- 
tion of properties, giving a constant contribution to A^eff 
but also late time energy density and acceleration. The 
interplay between these properties will allow each model 
to impact the observations in a distinctive manner. 

In addition to approaching this model microphysically, 
through the class of barotropic models, one can obtain an 
equivalent result within k-essence [33j using the quadratic 
Lagrangian £ = Aq + cA^, where A is the kinetic energy 
and c, Ao arc constants. 

Because for this model the relativistic scaling occurs 
so quickly (by ^ > 5), the AA^^^ contribution in this 
case will be constant at and before recombination, like 
the neutrino model with AAjf = 7.44 17f /(I - 17f ) (see 
Eq.[n]). However it has the late time change in equation of 
state that will affect large scale aspects of the CMB, and 
other cosmological probes, like the EDE model. Thus we 
expect the results to have aspects of each of the other 
two cases. 



D. Effects of the new components 



Since our main observational datasets are CMB 
anisotropics, it is useful to see the effects of sterile neu- 
trinos, EDE and barotropic dark energy on the CMB 
anisotropy angular spectrum. In Figure [2] (top panel) we 
show the CMB temperature angular spectra for these 3 
components assuming that they contribute at the level 
of one extra degree of freedom at BBN. As we can see, 
while sterile neutrinos and barotropic dark energy pro- 
duce nearly identical angular spectra, EDE predicts a 
significantly different spectrum. This is clearly shown 
in the bottom panel of the same figure where we plot 
the isolated Integrated Sachs Wolfe (ISW) contribution 
for each case. As we can see, the main difference be- 
tween a sterile neutrino and EDE comes from the ISW 
effect. This is mainly due to the time-dependence of the 
equation of state in the EDE component that tracks the 
dominant component at all epochs and increases the ISW 
signal on all angular scales. At the same time, we see that 
barotropic dark energy differs from a sterile neutrino in 
the increase in the ISW at large angular scales, due to 
the variation in the equation of state at small redshift in 
the barotropic component. 
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FIG. 2: CMB temperature (top panel) and ISW contribution 
alone (bottom panel) angular power spectra dependence from 
sterile neutrinos, early dark energy and barotropic dark en- 
ergy. All the models have been chosen to produce one extra 
relativistic degree of freedom at the epoch of BBN, except for 
the solid curve showing the standard case. 
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III. ANALYSIS AND RESULTS 

We perform a COSMOMC [s^] analysis combining 
the following CMB datasets: WMAP7 ACBAR |1, 
QUAD (3 (collectively referred to as "All"), and ACT Q. 
We analyze datasets using out to /max = 2500. We also 
include information on dark matter clustering from the 
galaxy power spectrum extracted from the SDSS-DR7 
luminous red galaxy sample [ssj . Finally, we impose 
a prior on the Hubble parameter based on the Hubble 
Space Telescope observations [36j |. 

The analysis method we adopt is based on the publicly 
available Markov Chain Monte Carlo package cosmomc 
[stI ] with a convergence diagnostic done through the Gcl- 
man and Rubin statistic. We sample the following nine- 
dimensional set of cosmological parameters, adopting flat 
priors on them: the baryon and cold dark matter densi- 
ties cjb and Wc, the Hubble constant Hq, the scalar spec- 
tral index ns, the overall normalization of the spectrum 
A a.t k ~ 0.05 Mpc~^, the SZ amplitude Asz, the optical 
depth to reionization, t, the effective number of relativis- 
tic neutrinos N^ff^, and finally the early density 51o, for ci- 
ther the case of EDE or the barotropic model (51^). For 
the ACT dataset we also consider two extra parameters 
accounting for the Poisson and clustering point sources 
foregrounds components. We consider purely adiabatic 
initial conditions and we impose spatial flatness. 

To study the impact of EDE perturbations, we con- 
sider two cases: "quintessence" (c^ = l,c^is = 0) and 
"relativistic" (c^ = 1/3 = c^jj EDE scenarios. For the 
barotropic dark energy model we assume = 1/3 = c^jg. 



A. Constraints on EDE with N^ff = 3.046 fixed 

We first perform an analysis of current data fixing the 
effective number of relativistic neutrinos to the standard 
value of Nl^ff = 3.046 and varying the amount of early 
dark energy, parametrized as He- The EDE affects the 
CMB angular power spectrum at all multipoles as shown 
in Fig. [5] We convert the EDE into an equivalent addi- 
tional relativistic species AN^^^ and quote this param- 
eter at the BBN epoch. We also recognize the fact that 
the changed expansion rate during the BBN, due to the 
presence of EDE, also alters the primordial Helium mass 
fraction Yp, and show the effect of ilg on Yp below. 

The MCMC results on the cosmological parameters are 
reported in the first two columns of Tabic H] 

As we can see, the cosmological data we consider do 
not provide evidence for an EDE component and signifi- 
cantly improve the bound |25| . yielding a 95% c.l. upper 
limit of rie < 0.043 in case of a "relativistic" EDE with 
Cg = cJig = 1/3, and a bound of il^ < 0.024 in case of 
a "quintessence" EDE with = 1, c^jg = 0. The as- 
sumptions of Cs and Cvis strongly affect the bounds on 
rig. The "quintessence" scenario leaves a stronger signal 
on the CMB anisotropics (see e.g. [l^]) and is therefore 
better constrained. This arises from the greater decay of 



potentials at recombination in contrast to the low sound 
speed case where the dark energy perturbations help sus- 
tain the potentials. 

In order to investigate the impact of the recent ACT 
dataset, which samples very small angular scales, giving 
a long lever arm, on the final result we also perform an 
analysis excluding it in the case of "relativistic" EDE. 
Without ACT we get a - 20% weaker bound, < 0.051 
at 95% C.L. 

While the EDE component is not preferred, it is also 
not excluded from current data. It is therefore interest- 
ing to investigate if the EDE component compatible with 
cosmological data is also compatible with BBN data. For 
this reason we translated the bounds on into the cor- 
responding bound on AN^^^ expected at time of onset 
of BBN, using Eq. ([5]), and computed the expected abun- 
dance Yp in primordial ^He by making use of the public 
available PArthENoPE BBN code (sec In other 

words, the constraints on early dark energy during the 
BBN correspond to limits on the expansion rate of the 
universe at this epoch, which translate into the corre- 
sponding limits on the excess of primordial mass fraction 
of Helium. 

Tabic U shows the "relativistic" case 95% upper limit 
f^e < 0.043 translates to a 95% constraint oi Yp ~ 
0.2504 ± 0.0026 (with ANf^f^ < 0.34), while the 
"quintessence" case 95% upper limit il^ < 0.024 trans- 
lates to a 95% constraint of Yp = 0.2495 ± 0.0016 (with 
AN^ff^ < 0.18). These values should be compared with 
the theoretical value of Yp = 0.2487 ± 0.0002 obtained 
assuming standard BBN and il,. = 0. EDE is therefore 
clearly shifting the BBN predictions on Yp towards larger 
values with weaker constraints. The weaker constraints 
indicate a degeneracy between f2e and Yp, as we discuss 
more in the next section, that we also show in Figure [3] 
where we plot the 68% and 95% constraints on the Yp- 
plane in the case of "relativistic" or a "quintessence" 
EDE. 

As stated in the introduction, current experimental 
measurements seems to prefer a larger value for the pri- 
or 



lilo: 
m [13 



mordial Hehum with Yp = 0.2561 ± 0.0108 (see 
Yp = 0.2565 ± 0.0010 (stat.) ±0.0050 (syst.) from 
These results are off by ~ 1.5(t from the expectations 
of standard BBN but introducing EDE acts to alleviate 
this tension. Given the possibility of systematics in mea- 
suring the primordial nuclear abundances, however, it is 
premature to derive any conclusion. 



B. Constraints on EDE and N^g 

As a second step, we include into the analysis the pos- 
sibility of extra sterile neutrinos, parametrizing it with 
the effective neutrino number N^f^. As we can see from 
Table [I] (last two columns), the constraints on flc are 
practically unaffected by the inclusion of extra RDOF 
and vice versa. From our analysis we found that sterile 
neutrinos are preferred with N^^g = 4.37 ± 0.76 at 68% 
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All+ACT 


Model: 


N^s = 3.046 


^^off varying 




Cs — ^vis — 1/3 


— 1 c^- — 


Cs — ^vis — 1/3 


— 1 c^- — 

^ y ^VIS ^ 


Parameter 










Ho 

Us 

to/ Gyr 

^cff 
fie 

AiVe^^^(aBBN) 


0.02218 ±0.00044 
0.1178 ±0.0039 
68.2 ± 1.7 
0.971 ±0.013 
13.71 ±0.30 
3.046 

< 0.043 

< 0.34 
0.2504 ±0.0013 


0.02232 ± 0.00044 
0.1163 ±0.0038 
67.8 ±1.6 
0.964 ±0.011 
13.83 ±0.29 
3.046 

< 0.024 

< 0.18 
0.2495 ± 0.0008 


0.02238 ± 0.00047 
0.138 ±0.012 

72.5 ±2.8 
0.988 ±0.015 
12.91 ±0.48 
4.37 ±0.76 

< 0.039 

< 0.32 
0.2661 ± 0.0078 


0.02259 ± 0.00048 
0.139 ±0.011 

72.4 ±2.7 
0.986 ±0.015 
12.94 ± 0.48 
4.49 ± 0.72 

< 0.020 

< 0.18 
0.2667 ±0.0080 



TABLE I: Best-fit values and 68% confidence errors on cosmological parameters using the current cosmological data. For 
and AN^^^(aBBN), EDE density and the contribution to the RDOF from EDE at the BBN epoch respectively, the upper 
bounds at 95% c.l. are reported. See text for other details. 
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FIG. 3: 68% and 95% c.l. contours in the Yp-Qc plane for 
the standard EDE model. The red dashed contours show 
the Cg = Cvis = 1/3 model, while the blue solid contours 
show the Cs = l,Cvis = model. Since the early dark energy 
enhances the expansion rate during the BBN, it allows for a 
higher primordial Helium mass fraction according to Al^ ~ 
0.013(iVc£f - 3) [3. 



EDE behaves more like matter and so can be constrained 
separately from the neutrino contribution. 

This can also be seen in Figure |4] where we show 
the likelihood contour plots in the N^g-ftc plane for the 
cases of "relativistic" and "quintessence" EDE. There is 
no strong degeneracy between J7e and N^g. This, to- 
gether with the small value of AN^^^ allowed, indi- 
cates that the current hints for the existence of the extra 
RDOF cannot be completely explained by a conventional 
EDE model. In the next subsection we will see that the 
barotropic class of early dark energy has more success. 

As we can see from Table HI including the possibil- 
ity of extra neutrino contributions to TVoff greatly en- 
larges the CMB bounds on primordial Helium, with 
Yp = 0.2661 ± 0.0078 in case of "relativistic" EDE to 
Yp = 0.2667 ± 0.0080 in case of "quintessence" EDE. 
This stronger influence of neutrino RDOF than EDE has 
the consequence that in this situation the impact of EDE 
on the Yp abundance is small. As seen from the results in 
Table m the AN^ff^ from EDE at BBN is always better 
constrained from CMB data than the Nl^g expected from 
a sterile neutrino. If future measurements of primordial 
^He clearly point towards value of Yp ~ 0.26, it will not 
be possible to explain this result with a conventional EDE 
contribution. 



C.L. This constraint should be compared with the bound 
from the analysis of ^ of N^g = 4.6 ± 0.8 at 68% c.L, 
obtained with similar datasets but without EDE, indi- 
cating that the effect of EDE on the constraint is small. 
The low covariance between the number of sterile neutri- 
nos and EDE density comes from the property that while 
at BBN they both act as RDOF, by recombination the 



Finally, we note that including the possibility of N^^^ > 
3 also changes the constraints on making it more com- 
patible with a Harrison-Zeldovich, = 1, primordial 
spectrum (cf. [13, H^)- The best fit value of the infiation- 
ary tilt ris — 1 is reduced by almost a factor of 3, which 
would have a substantial impact in the reconstruction of 
the inflationary potential. 
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Parameter 


All + ACT 




0.02209 ± 0.00055 




0.135 ±0.012 


Ho 


71.1 ± 2.8 


Us 


0.986 ±0.015 


to/ Gyr 


13.18 ±0.51 




< 5.1 




< 0.37 




< 2.8 




0.2649 ± 0.0084 



0.05 



TABLE II: Best-fit values, together with 68% confidence er- 
rors, on cosmological parameters for the barotropic model us- 
ing current data. For the A'efi, and AA^^ parameters 
upper bound at 95% c.l. are reported. 



FIG. 4: 68% and 95% c.l. contours in the N^s-D.c plane for 
the standard EDE model plus neutrinos). The red dashed 
contours refer to Cg = 0^;^ = 1/3 case, while the blue solid 



contours refer to the c. 



case. 



C. Results on Barotropic Dark Energy 

The barotropic model contributes both early dark en- 
ergy density and a constant (rather than diluted as in 
Fig. [T]) early time RDOF. This will have interesting im- 
plications. As for the conventional EDE case, we add 
the early density, here il^ , to the MCMC analysis to es- 
timate constraints on cosmological parameters. We also 
allow N^ff to vary as in the previous section. For sim- 



plicity, we otherwise fix wq = —0.99 and 



1/3. 



The results are reported in Table [III and in Figure [S] we 
show the degeneracy between TV^g and flf parameters. 

The barotropic model strongly alters the constraints on 
N^ff and a non-negligible presence of the dark radiation 
part of the barotropic dark energy at recombination could 
not only bring the constraints on N^g back in agreement 
with the standard value of A^^fj = 3.046 but even erase the 
current claim for a neutrino background from CMB data. 
A "neutrinoless" model with N^^^ = and il^ = 0.4, 
albeit extreme, is allowed by the cosmological data we 
consider here. 

As in the case for N^^^f > 3, when a barotropic dark 
energy model is considered (even without extra neutri- 
nos) a high value of Yp is consistent and the constraints 
on rig are moved toward a Harrison-Zeldovich primordial 
spectrum. 



D. Forecasts for the Planck Satellite Mission 




0.5 



FIG. 5: 68% and 95% c.l. contours in the N^g-^c plane for 
the barotropic dark energy model. 



achievable from this satellite, we consider an experimen- 
tal configuration with three frequency channels with the 
specifications as listed in Table IIIII (see [i^ ) . 



Experiment ChannelfGHz] 


FWHM 


ar [fJ-K] 


C7p[llK] 


Planck 143 


7.1' 


6.0 


11.4 


/sky = 0.85 100 


10.0' 


6.8 


10.9 


70 


14.0' 


12.8 


18.3 



TABLE III: Planck experimental specifications. 



Looking to the future, we investigate the constraints on 
EDE and RDOF in the case of the data as expected from 
the Planck satellite. To evaluate the future constraints 



For each frequency channel we consider a detector noise 
of (6*0-)^ where 6 is the FWHM of the beam assuming a 
Gaussian profile and a is the sensitivity. We therefore 
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take a noise spectrum given by 

7V/ = (0^^)2ei(;+i)/i? 



(10) 



where lb = VS In 2/6' and the label X refers to either 
temperature or polarization, X = T, P. 

We perform the standard Fisher matrix analysis eval- 
uating (see e.g. (HI): 



(11) 



where £(data|p) is the likelihood function of a set of pa- 
rameters p given some data, and the partial derivatives 
and the averaging are evaluated using the fiducial val- 
ues po of the parameters. The Cramer-Rao inequality 
implies that {F~^)ii is the smallest variance in the pa- 
rameter Pi, so we can generally think of F^^ as the best 
possible covariance matrix for estimates of the vector p. 
The one sigma error forecasted for each parameter is then 
given by cTp^ = y^{F-^)u. 

We consider a set of 10 cosmological parameters 
with the following fiducial values: the physical bary- 
onic and cold dark matter densities relative to criti- 
cal ribh'^ = 0.02258 and flch'^ = 0.1109, the optical 
depth to reionization r = 0.088, the Hubble parameter 
Hq = 71km/s/ Mpc, the current dark energy equation 
of state Wo = —0.95, the early dark energy density rela- 
tive to critical fio = 0.03, the spectral index Us = 0.963, 
and the number of relativistic neutrinos N(.e = 3.046. 
For the last two parameters, the effective and viscous 
sound speeds and c^jg, we choose alternate fiducial val- 
ues of (1/3,1/3) (the "relativistic" model) or (1,0) (the 
"quintessence" model). 

1. EDE Forecasts 

In Table IIVI we report the uncertainties obtained 
on the cosmological parameters. The degeneracy 
between fig and N^^^ is shown in Figure |6] for the 
two analysed cases. As seen in the Figure and in 
the Table, the future data from Planck will provide 
strong constraints on the RDOF: cr(A'cff) — 0.11, 
with little impact from the EDE density. If EDE 
with fie = 0.03 is present, it will be detected at high 
significance, since (T(rie) ~ 0.005. The radiation 
and quintessence configurations of EDE can also 
be distinguished. 

2. Barotropic DE Forecasts 

Similarly to the previous analysis, we forecasted 
the errors on cosmological parameters with data 
expected from Planck in a barotropic dark energy 
scenario. We report in Table IVl the l-cr errors, and 
in Figure [7] we show the degeneracy between N^^^ 
and nf. 

For a fiducial early density of = 0.03, the 
barotropic model cannot be readily distinguished 
from the standard EDE model. However, if ftf^ is 





Planck 1-a 


uncertainty 


Model: 


~ Cvis ~ 1/^ 


Cg 1 , Cyis 


Parameter 


Fiducial 








0.02258 


0.00016 


0.00014 




0.1109 


0.0018 


0.0017 


T 


0.0880 


0.0020 


0.0022 


Ho 


71.0 


8.5 


8.8 


ris 


0.9630 


0.0046 


0.0044 




3.046 


0.11 


0.11 


wo 


-0.95 


0.24 


0.24 


fie 


0.030 


0.005 


0.004 




0.33 


0.047 




'-'VIS 


0.33 


0.13 






1.00 




0.34 


^vis 







0.11 



TABLE IV: Fiducial errors and forecasted 1-a errors expected 
from the Planck satellite in the EDE scenario. 



3.50 
3.25 

3.00 

2.75 
2.50 




0.01 0.02 0.03 0.04 0.05 

0. 



FIG. 6: 68% and 95% c.L constraints in the N^g-Qe plane. 
Solid lines show the "relativistic" case = c' 



1/3, while 
= 0. 



dashed lines show the "quintessence" case — 1, 
The fiducial values are given by the "+" symbol. 



much larger, then distinction will be possible, with 
the associated implications for RDOF, Yp, and n^. 



IV. CONCLUSIONS 

Data analyses of recent cosmological data have re- 
ported an interesting indication for the presence of an 
extra background of relativistic particles. In this paper 
we have investigated the stability of this result by consid- 
ering the influence of a possible early dark energy com- 
ponent. We have found that current data do not provide 
evidence for an EDE component, updating and strength- 
ening previous constraints on EDE, although there is still 
room for an interesting contribution. In particular, we 
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Planck 1 — cr uncertainty 


Model: 


Cs ~ Cvis ~ 


Parameter 


Fiducial 






0.02258 


0.00013 




0.1109 


0.0019 


T 


0.0880 


0.0022 


Ho 


71.00 


0.88 


ris 


0.9630 


0.0041 




3.046 


0.17 


Wo 


-0.95 


0.041 




0.030 


0.015 




0.33 


0.045 


Wis 


0.33 


0.17 



TABLE V: Same as Table |IV1 but for the barotropic dark 
energy scenario. 




0.000 0.025 0.050 0.075 0.100 



FIG. 7: 2-D contour plots at 68% and 95% c.l. in the plane 
N^^ff-Qf for Cg — Cvis = 1/3. The fiducial values are reported 
with the "+" symbol. 



acy between them and the constraints are virtually unaf- 
fected. The indication for extra neutrinos in current data 
at about two standard deviations is unchanged even after 
allowing a EDE component. We found Nj^^f = 4.37 ±0.75 
at 68% c.l. for relativistic EDE and N^^^ = 4.49 ± 0.72 
for quintessence EDE when CMB, SDSS-DR7 and HST 
data are combined. The bounds on are practically 
unchanged. The key point is that EDE starts to behave 
differently from a relativistic component after radiation- 
matter equality, before the epoch of recombination. CMB 
data can therefore provide crucial information in discrim- 
inating between N^g and early dark energy while for BEN 
these two components are virtually indistinguishable. 

However, when a barotropic dark energy model is 
considered, we have found that the constraints on Nj^g 
can be strongly altered, bringing the standard value of 
Nl^ff ~ 3.046 back into perfect agreement with observa- 
tions. In fact, even a "neutrinoless" model with N^g = 
and = 0.4 is allowed given current observations. 
While that is extreme, the model dependency clearly in- 
dicates the caveats of considering A^cft > 3.046 as an indi- 
cation for an extra sterile neutrino or claiming any detec- 
tion for a neutrino background. Interestingly, barotropic 
dark energy also shifts Yp to higher values, and reduces 
jris — 1|, with implications for models of inflation. 

Finally, we have shown that for the Planck experi- 
ment alone, again no substantial degeneracy is expected 
between f2e and N^g, with an expected accuracy of 
a{N^g) = 0.11 and ai^e) ^ 0.005. However, in a 
barotropic dark energy scenario the degeneracy is present 
and Planck will offer weaker bounds with an estimated 
a{N^g) = 0.17 and cr(f2f ) = 0.015. Both are still a great 
improvement over present data. Planck will therefore 
shed light on various scenarios for early dark energy and 
relativistic degrees of freedom, exploring if new physics 
exists in the neutrino or dark energy sectors, a possible 
shift in the inflationary tilt rig — 1, and the consistency 
of the primordial Helium abundance Yp. 



found the following 95% c.l. upper limits: fie < 0.043 for 
"relativistic" EDE (c^ = c^;, = 1/3) and fie < 0.024 for 
"quintessence" EDE (c^ = 1, c^jg = 0). These bounds 
translate into an extra relativistic background at BBN of 
AN^g^^ < 0.34 and ANgf"^ < 0.18 at 95% c.l., respec- 
tively. 

The EDE models are therefore not able to change the 
amount of primordial ^He produced in BBN by more than 
AYp^^ = 0.0044, and do not help much in explaining 
why recent measurements of abundances of primordial 
Helium show values larger than those expected from stan- 
dard BBN. The systematics in those measurements are 
however still too large to conclude that there is conflict 
between the measured and predicted "^He abundance. 

When both an EDE and extra sterile neutrinos are 
considered in the analysis, there is very little degener- 
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